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SUMMARY 


A constant-property analysis of the heat transfer at a nonisothermal surface of a blunt 
body has been made by means of series expansions of Blasius type for two dimensional and 
axially symmetrical laminar flow. The temperature of the surface and the velocity at the 
edge of the boundary layer are assumed to have power series distributions. Some new terms, 
containing universal functions, must be added to the existing solutions for heat transfer at 
isothermal surfaces, earlier published by the author. Some of these functions have been cal- 


culated numerically. 


1. Introduction 


The problem of calculating the heat transferred by forced convection from 
a body with a temperature varying over its surface, can arise under several 
circumstances. Thus, even at stationary conditions, a nonuniform heating or 
cooling from inside or outside of the body can result in a nonuniform tem- 
perature distribution owing to the limited heat conductivity of the matter. 
However, whether the real temperature of the surface becomes practically uni- 
form or not, the problem mentioned will arise if the velocity of the flow is so 
large that heating effects due to internal friction and to pressure changes are 
noticeable. At not too large velocities the energy equation describing the tem- 
perature field remains linear, the inhomogeneous terms giving account of those 
heating effects. A general solution can then be obtained by adding a particular 
solution of the complete equation and a general solution of the homogeneous 
part. This homogeneous part is the equation of low speed heat transfer [1]. 
For the former solution, that of the adiabatic wall (that is, that of a non- 
heated, insulated body) is useful. The total heat transfer will then be described 
by the homogeneous equation and will depend on the difference between the 
actual wall temperature and that of the adiabatic case. At moderate velocities 
the properties of the fluid can be treated as constants [2]. 

Much work has been done on exact solutions of the heat transfer at non- 
isothermal surfaces of plates, wedges and cones (e.g. [3], [4], [5]). For cal- 
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culating boundary layers at bodies of arbitrary shape, there are available several 
approximate methods which, however, require a further assumption together 
with the fundamental equations. For getting solutions without such assumptions, 
a series expansion with respect to the distance along the wall can be made. 
It can be continued by some step-by-step method. For the thermal boundary 
layer at vanishing viscous heating and isothermal surfaces, such a series solu- 
tion has been worked out for blunt bodies by the author [6] for two-dimensional 
and axially symmetrical flows. In the former case the solutions given by Blasius [7], 
Hiemenz [8], and Howarth [9] were used for the velocity field. The corre- 
sponding solutions for the viscous heating at adiabatic walls have been studied 
by the author [10]. The series solution given by [6] can be extended to a 
constant-property solution in the nonisothermal case, still using universal func- 
tions independent of the particular data of the problem. A brief account of 
this is given below. The transition from the dimensioned quantities to the 
dimensionless universal functions has now been made in a way somewhat dif- 
ferent from that of paper [6]. Since the equations of heat and mass transfer 
are quite the same in this case, the analysis can be used also for mass transfer. 
Stationary conditions are assumed. 


2. Series expansions 


A. Two-dimensional symmetrical case 


The velocity- and temperature fields are described by the equations of mo- 
mentum, continuity, and energy; and by the corresponding boundary and initial 
conditions. For a laminar boundary layer the energy equation for a constant- 
property analysis at neglected viscous heating is 

U ing Obit, Cat é 
Ou dy dy?’ (1) 
(x, y and wu, v=coordinates and velocity components parallel and perpendicular 


to the wall; a= thermal diffusivity =A/oc,; t=temperature excess with respect 
to the surroundings). The velocity U outside the boundary layer is assumed to be 


U =u, 0+ Ugx? +ugx? ++ (2) 
(the leading term indicating a blunt body). The boundary conditions are: 
if y=0, t prescribed; and if y= oo, t=0. (3) 


For the isothermal case (prescribed value =constant =t)) ref. [6] gives the 
solution 


a 2 4 
Cot Cyt + Cy 2+ ++ (4) 
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The quantities with capital letters are universal functions depending (besides 
on the Prandtl number o=v/a) only on n=yVu,/r; (v=kinematic viscosity). 
For those functions the following equations are valid. 


taal § ; 

- Fo = -f,Fo 
1 47 , , / 

att =F s + 2 fF, = 35 Fo 

1 y// j / , 

mp = —f,G4+4fiG,—595 Fo 
1 ‘r , , y 8 , / | 
Bot = ~f Het 4h Hy— dhs Pot 3 (2fsP,—3f, Fs) t (6) 
it yft / , / | 
578 = —/,G,1+6f/1:6,—79,Fo | 
1 v7 , / , , , / / 
ais = —f,He+6fHe—Th, Fo+3(4fsG,+ 29; F,—3f,G1—595 Fo) ‘ 


1 rir , , , / , , , 
pHs aa —f, Ket 6fi:Kk,—7k, Fo+3(4feH,+2hs F,—3f,Hs— 5h; Fo) 


JJ 


The functions f,... (small letters) are the velocity functions of Howarth [9]. 
The boundary conditions of the equations (6) are 7=0, F,=1 and other func- 
tions=0; 7= co, all functions=0. Some of the functions are calculated numeri- 
cally and given in [6]. 

The heat transfer is described by the Nusselt number Nu. If D and t,, are 
the reference quantities for getting a dimensionless number, we have 


: (7) 


Nu can be calculated from equations (5), (6), (7). 
The solution described can be generalised to the nonisothermal case, if the 
prescribed temperature distribution over the surface is of power series type: 


t=t, tigv* tia*+ ++. (8) 


Of course the equations (6) could be solved directly by satisfying the par- 
ticular coefficient values of (8), but this would not give any universal functions 
independent of the particular data of the problem. Instead of using this way 
we benefit by the fact that the energy equation now is linear and therefore 
solutions can be superimposed. We therefore study for »>0 the distribution 


t=t,2" (nm an even number). (9) 


Since the equations of the velocity field are not linear we cannot superimpose 
solutions made for the single terms of (2), that is for wedge solutions. Now 
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for studying the distribution (9) we can still use (5) and (6), together satisfying 
(1), only changing in (5) the constant ¢, for 


“4 (10) 


t 
n 
Un+1 


and giving some of the functions new boundary values at the wall. Thus for 7=0, 
the first function in the factor for x” gets the value 1/(m+ 2) and all other func- 
tions vanish. The former function is called C,,. Some of the latter ones are defined 
by homogeneous equations and will therefore, according to the boundary values, 
become identically equal to 0. Those that do not disappear but possibly get 
their equation slightly changed, are marked by a circle (for example K,). 

By summing solutions for different n-values we get the solution for the distri- 
bution (8): 


4u ° 
=, Poses ye? eee +t, |40,027+6 2H, xt + 
Uy ; ? 068. ine 
: ut 


2 
U ° U ° ° 
+8 (“H+ | a + -.| +t, locas =3 Fy 2° + -.| . 
Uy Uji U 


1 


tg [8 Cyah bm |e (11) 


Of course the terms in (11) could be collected according to powers of x. The 
functions are defined by the equations 


1 ut , , 
me = —f,C2+2fiC, 


1 ” , , 
get TA Oat 4h, 
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| 
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H=oco, all functions=0. 
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The final expression for Nw is obtained as follows: 


D Uy , 4u , , ot 
Nu =-— ma |Fo+ 88 Rats + te 4030? + 6-9 Hiat + 
tm v Uy oi 


2 
Un °, Ug 2, , u ca 
+8(“8 ri. + 43 he) a + | i Js o, 24 + 82 Fi a + -.| - 
Uy Uy Uy 


+t, [BCpaP + ]4 | : (13) 


n=0 


The set of terms containing ¢, is the solution for the isothermal case, as given 
in [6] by the author. 


B. Axially symmetrical case 


Owing to the different form of the equation of continuity in this case some- 
what different formulas are obtained. The radial distance r is supposed to be 
of the form 


r= C+ ee +75 + (14) 


and the velocity distribution of the form (2). 
Using the equations of ref. [6] we get for a wall temperature distribution 
of the type (8) in analogy with the case A the functions defined by the fol- 


lowing equations (now 7=y V2u1/y): 


1 tr , ) 4 

i: = —f,Cot+fic, 

1 7 / , 

~ C6 = —f,,+2fiC, | 

files dy met Ose hee PAS | 

o 3 1 4 a 2 | (15) 
1 ony eh, ae 2 , , , 

pet = Thidet 2h Tat 3 (2hsCy— Aha Ca— fr Ca) | 


Toy , D 
58 = —f,Cet+3fiC, 


2 
n=0, C,=—; and other functions =0. 
Boundary conditions: n+2 


n=, all functions=0. 
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Final expression for Nw: 


D Qu , 2u , rg U / 
Nu=—-—— /— I, [i+ 288(as +f y;) eet 
bin v \ Uy 1, Us 
3u ru UR rau ru 
ae (ere Ca S Kit 3 3h4+5 1 Qi) ate | + 
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+t [4Cga°+ ---]+ lye (16) 


3. Numerical calculations 


In the two-dimensional case numerical calculations were made for the seven 
universal functions of (12) with a Prandtl number o=0.7 (air). 

In the axially symmetrical case the corresponding calculations were made for 
the first four and five functions of (15) with the Prandtl numbers o=0.7 and 
1/0.395 respectively (the latter is equal to the Schmidt number for naphthalene 
and is thus applicable to the evaporation of this substance). 

For the velocity functions the same values as in [6] were used. Some groups 
of the equations had common homogeneous parts which were equal to those 
of certain groups of the equations in [6] and [10]. This made simplifications 
possible. For such a homogeneous equation one particular solution was computed 
by the Runge-Kutta method. A second particular solution, linearly independent 
of the former, was obtained with the aid of the Wronski determinant, which 
gave a first order differential equation for this second solution. This was possible 
since the Wronski determinant can be deduced from the next to the highest 
term [11]. The final solution of the homogeneous equation was obtained by 
linear combination of those two particular ones. A particular solution of an 
inhomogeneous equation in the group could be computed by variation of the 
constants, which could be done by numerical quadrature. Thus an equation of 
the type 


z’+of,2’-—mofiz=f (17) 


(f, and f functions of 7) with a particular solution z, for the homogeneous part 
and with the actual boundary conditions has the solution 


oo t 
Ff 
aa Pf [a(| Fes)ac 2 
ema | 5({ Fae)ae—a' - | Sac: (18) 
0 0 | gas 
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where F=e 9° (19) 


For the particular solution z, it was convenient to use the initial conditions 
n=0; z,=arbitrary value, z;=0. 

The results obtained for the functions are shown in Tables 1 and 2 respec- 
tively and in Figs. 1 and 2. The final decimal place is correct for C,, and has 
an error of at most a few units for the other functions. 


The following values of the derivatives were obtained at the wall (7 =0) with 
all decimal places correct: 


Two-dimensional : 
o=0.7; 0;=—0.2130; Cy=—0.1757; Hi=—0.1113; C= —0.1505; 
Hoa= — 0.1096: Hi, = —0.1006; Ki = + 0.0414, 


Axially symmetric: 
o=0.7; C= — 0.3333; C= —0.2657; Ji=—0.0323; Ki= —0.1653; 
o=1/0.395; Cz= — 0.5306; C,= —0.4203; Ji= —0.0497;  Ki= — 0.2787. 
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Table 1. Thermal universal functions for two-dimensional symmetrical case. 6 = 0.7. 


n C; CO, HA, C, Hea ' Heo 1G 


0 0.2500 0.1667 0 0.1250 0 0 0 
0.2 0.2079 0.1322 — 0.0214 0.0956 — 0.0211 — 0.0190 0.0084 
0.4 0.1685 0.1010 — 0.0391 0.0698 — 0.0385 — 0.0329 0.0171 
0.6 0.1334 0.0748 — 0.0511 0.0490 — 0.0502 — 0.0402 0.0256 
0.8 0.1033 0.0538 — 0.0574 0.0333 — 0.0562 — 0.0417 0.0332 
1.0 0.0784 0.0378 — 0.0586 0.0220 — 0.0571 — 0.0390 0.0390 
1.2 0.0583 0.0260 — 0.0559 0.0142 — 0.0542 — 0.0340 0.0428 
1.4 0.0425 0.0174 — 0.0506 0.0089 — 0.0487 — 0.0281 0.0444 
1.6 0.0304 0.0115 — 0.0438 0.0055 — 0.0420 — 0.0221 0.0437 
1.8 0.0213 0.0074 — 0.0365 0.0033 — 0.0347 — 0.0169 0.0415 
2.0 0.0147 0.0047 — 0.0294 0.0020 — 0.0279 — 0.0123 0.0376 
2.2 0.0099 0.0029 — 0.0230 0.0012 — 0.0215 — 0.0089 0.033 
2.4 0.0065 0.0018 — 0.0174 0.0006 — 0.0163 — 0.0064 0.028 
2.6 0.0042 0.0011 — 0.0128 0.0004 — 0.0118 — 0.0043 0.023 
2.8 0.0027 0.0006 — 0.0091 0.0002 — 0.0085 — 0.0027 0.018 
3.0 0.0016 0.0004 — 0.0064 0.0001 — 0.0061 — 0.0017 0.016 
3.2 0.0010 0.0002 — 0.0043 0.0000 — 0.0041 — 0.0011 0.011 
3.4 0.0006 0.0001 — 0.0028 — 0.0026 — 0.0008 0.008 
3.6 0.0003 0.0000 — 0.0018 — 0.0016 — 0.0005 0.005 
3.8 0.0002 — 0.0011 — 0.0009 — 0.0003 0.004 
4.0 0.0001 — 0.0007 — 0.0004 — 0.0001 0.002 
4.2 0.0000 — 0.0004 — 0.0002 — 0.0000 0.001 
4.4 — 0.0002 — 0.0000 0.001 
4.6 — 0.0001 0.000 
4.8 — 0.0001 

5.0 — 0.0000 


Table 2. Thermal universal functions for axially symmetrical case. o = 0.7. 


ie Nea, | C, | di ke n | C, | C, | Te | K, 

|| 
0 | 0.5000 | 0.3333 0 0 | 3.2 | 0.0054 | 0.0017 | —0.0066 | — 0.0188 
0.2 | 0.4338 | 0.2807 | —0.0064 | —0.0325 | 3.4 | 0.0034 | 0.0010 | —0.0047 | —0.0132 
0.4 | 0.3698 | 0.2309 | —0.0126 | —0.0617 | 3.6 | 0.0021 | 0.0006 | —0.0033 | —0.0090 
0.6 | 0.3098 | 0.1857 | —0.0183 | —0.0855 | 3.8 | 0.0013 | 0.0003 | —0.0022 | —0.0060 
0.8 | 0.2550 | 0.1463 | —0.0230 | —0.1027 | 4.0 | 0.0007 | 0.0002 | —0,0014 | —0.0038 
1.0 | 0.2062 | 0.1129 | —0.0266 | —0.1128 | 4.2 | 0.0004 | 0.0001 | —0.0009 | —0.0024 
1.2 | 0.1637 | 0.0854 | —0.0288 | —0.1162 | 4.4 | 0.0002 | 0.0000 | —0.0006 | —0.0014 
1.4 | 0.1275 | 0.0633 | —0.0296 | —0.1137 | 4.6 | 0.0001 — 0.0003 | —0.0008 
1.6 | 0.0975 | 0.0460 | — 0.0290 | —0.1065 | 4.8 | 0.0000 ~ 0.0002 | —0,0005 
1.8 | 0.0731 | 0.0328 | —0.0272 | —0.0960 | 5.0 — 0.0001 | —0.0002 
2.0 | 0.0537 | 0.0229 | —0.0246 | —0.0836 | 5.2 —0.0000 | —0.0001 
2.2 | 0.0387 | 0.0157 | —0.0215 | —0.0703 || 5.4 — 0.0000 
2.4 | 0.0273 | 0.0105 | —0.0181 | —0.0573 
2.6 | 0.0188 | 0.0069 | —0.0147 | —0.0453 | 
2.8 | 0.0127 | 0.0045 | —0.0116 | —0.0348 | 
3.0 | 0.0084 | 0.0028 | —0.0089 | — 0.0259 | 
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